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Abstract 

We have performed a coupled channel calculation of the 1 ++ cc sector including qq and DD* 
molecular configurations. The calculation was done within a constituent quark model which suc- 
cessfully describes the meson spectrum, in particular the cc 1 sector. Two and four quark 
configurations are coupled using the 3 Po model. 

The elusive X(3872) meson appears as a new state with a high probability for the DD* molecular 
component. When the mass difference between neutral and charged states is included a large 
D°D*° component is found which dominates for large distances and breaks isospin symmetry in 
the physical state. The original cc(2 3 Pi) state acquires a sizable DD* component and can be 
identified with the X(3940). We study the B -> Ktt + tt~ J/ip and B -)■ KD°D*° decays finding a 
good agreement with Belle and BaBar experimental data. 
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I. INTRODUCTION. 



In the last years a number of exciting discoveries of new hadron states have challenged 
our description of the hadron spectroscopy. One of the most mysterious states is the well 
established X(3872). It was first discovered by the Belle Collaboration in the J/ipirK 
invariant mass spectrum of the decay B + — > K + -k + 'k~ J/ip Its existence was soon 
confirmed by BaBar [2], CDF 3] and DO {4] Collaborations. The world average mass 
is M x = 3871.2 ± 0.5 MeV and its width T x < 2.3 MeV. The measurements of the 
X(3872) 7 J/V> decay implies an even C-parity. Moreover angular correlation be- 

tween final state particles in the X(3872) — > n + n~ J/ip decay measured by Belle [5J] suggests 
that the J PC = ++ and J PC = H may be ruled out and strongly favors the J PC = 1 ++ 
quantum numbers although the 2 ++ combination cannot be excluded. A later analysis by 
CDF Collaboration |7j of the same decay is compatible with the Belle results and concludes 
from the dipion mass spectrum that the most likely quantum numbers should be J PC = 1 ++ 
but cannot totally exclude the J PC = 2~ + combination . These conclusions were confirmed 
by a new CDF analysis of the decay X(3872) — > 7i + 7i^J/ip followed by J hb — > /U + /i~ ex- 
cluding all the other possible quantum numbers at 99.7% confidence level [8j. However the 
small phase space available for the decay X(3872) -> D°D°n observed by Belle Q discards 
the J = 2 leaving the 1 ++ assignment as the most probable option. 

In the 1 ++ sector the only well established state in the PDG 10j is the Xci(l-P) with 
a mass M = 3510.66 ± 0.07 Me V. The first excitation is expected around 3950 MeV. In 
this energy region Belle has reported the observation of three resonant structures denoted 
by X(3940_K F(3940) and Z(3930). The last one was observed by Belle in the 77 DD 
reaction [ll| and is already included in the PDG as the Xc 2 (2-P). The X(3940) has been 
seen as a peak in the recoiling mass spectrum of J/ip produced in e + e~ collision. Its main 



decay channel is DD" 



12J. The F(3940) appears as a threshold enhancement in the J/ipu 



invariant mass distribution of the B — > J/tjjuK decay 13] . 

The relative decay rates outlines a puzzling structure for the X(3872). The jJ/ip and 
'yip' decay rates fl^ 



X(3872) -» jj/ip 
X(3872) ii+Ti-J/i) 

X(3872) -» 7^ 
X(3872) Ti+ii-J/i) 



0.33 ±0.12 
1.1 ±0.4 



(1) 
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suggest a cc structure whereas the X(3872) — > n + n ir J/ip decay mode 

X(3872) ->■ n + TT-ir°J/ip 



X(3872) ir+<K-J/ip 



1.0 ±0.4 ±0.3 



(2) 



indicates a very different one 151 ]. The dipion mass spectrum in the 7r + 7r~ J ' /ip channel shows 
that the pions come from the p° resonance. On the other hand the 7r + 7r~7r° mass spectrum 
has a strong peak around 750 MeV suggesting that the process is dominated by a a; meson. 
Thus the ratio R ~ 1 indicates that there should be an isospin violation incompatible with 
a traditional charmonium assumption. 

Concerning the mass value, in 2006 Belle measured an enhancement in the D D°7i° 
channel just above the D°D*° threshold using the B + — > K + D°D°7r° decay. The amazing 
aspect of this enhancement is that it appears at M x = 3875.2 ± 0.7t?;| ± 0.8 MeV just 
3 MeV above the Mx world average mass value. This fact triggered a new discussion about 
the possibility of two different charmonium like states. The Belle mass value was confirmed 

n 

later by the BaBar Collaboration [16]. Last year the Belle Collaboration announced a new 



171 ] with a lower position of the X(3872) peak 



measurement of the B — > KD°D°tt° decay 
in M x = 3872.6±o;4 ± OA MeV. New data of the 7i + 7i~J/^ decay has been also recently 



reported by the Belle 



18|], BaBar [19| and CDF 



20| Collaborations, confirming a mass value 



in agreement with the world average. 



The X(3872) mass is difficult to reproduce by the standard quark models (see Ref. 21] 
for a review). The state appears to be too heavy for a ID charmonium state and too light 
for a 2P charmonium one. Moreover no four-quark bound state configurations have been 
found in this mass reg ion which rules out the possibility that this particle was a compact 
tetraquark system 22j, ]23|] . 

An important property of the X(3872) is that its mass is extremely close to the D°D*° 
threshold with a difference using PDG values given by — 0.6±0.6 MeV. The proximity of the 
D°D*° threshold made the X(3872) a natural candidate to a C = + D°D*° molecule. The 
hypothesis of a DD* molecule mainly bound by pion exchange has been suggested by several 
authors [24j. In particular, in Ref. 2M lt is argued that the X(3872) is a J PC = 1 ++ D°D*° 
molecule stabilized by admixture of pJ/ip and uJ/ip states. The author shows that pion 
exchange alone can not bind the molecule being the combined effect of pion exchange and 
coupled channels responsible for that. The D°D*° component dominates the wave function 
at the experimental binding becoming all other contributions small. 
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The molecular interpretation runs into trouble when it tries to explain the high jiff' 
decay rate. For a molecular state this can be only proceed through annihilation diagrams 
and hence is very small. 

This puzzling situation suggests for the X(3872) state a combination of a 2P cc state and 
a weakly-bound D°D*° molecule Q, Q]. The experimental assignment J PC = 1 ++ favors 
this conclusion because it allows the molecule to be in a relative S'-wave state whereas the 
corresponding cc should be in a relative P-wave state. Then the masses of the additional 
light quarks are compensated by the angular momentum excitation and both configurations 
may be almost in the same mass region. Similar behavior has been already observed in 
the open charmed sector [26|. Recently Zhang et al. [27| have analyzed, using the coupled 



channel Flatte formula, the B -»■ KD°D 7r° [l7fl and B -> K-k + -r- J/tf [l8j Belle data. They 
found that a third sheet pole close, but below, the D°D*° threshold is needed to describe the 
data, which supports the idea of the X(3872) as a mixed state of x' c i an d D°D*° components. 
An updated Flatte analysis of the same data together with the new BaBar data of the same 
reactions [3] has been performed in Ref. assuming a mechanism for the X(3872) 
production via the charmonium components. The authors conclude that the data clearly 
indicates a sizable cc 2 3 Pi component in the X(3872) wave function. Finally Dong et al. 29 ] 
show in their analysis of the J/ip 7 and i/j(2S) 7 decay modes of the X(3872) that the large 
value of the ratio PP(X(3872) -»■ J/^ 7 )/PP(X(3872) ip(2S)j) measured by the BaBar 
Collaboration provide a constraint on the value of the cc component in the X(3872). From 
the experimental values they deduce a small admixture of the cc component 

Having in mind these evidences, in this paper we perform a microscopic coupled channel 
calculation of the 1 ++ sector including both cc and DD* states. The calculation is done 
in the framework of a constituent quark model widely used in hadronic spectroscopy. The 
paper is organized as follows. In the next section we review the main ingredients of our 
model. Section III is devoted to discuss the numerical procedures and the results . Finally 
we summarize the main findings of our work in the last section. 
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II. THE MODEL. 



A. The constituent quark model 

The constituent quark model used in this work has been extensively described elsewhere 



30 



3 11 ] and therefore we will only summarize here its most relevant aspects. The model is 
based on the assumption that the light constituent quark mass appears as a consequence of 
the spontaneous breaking of the chiral symmetry at some momentum scale. As a consequence 
the quark propagator gets modified and quarks acquire a dynamical momentum dependent 
mass. The simplest Lagrangian must therefore contain chiral fields to compensate the mass 



term and can be expressed as [32 1 



£ = ip(ip- M{q 2 )W & )ij) (3) 

where U 15 = exp(iir a \ a j 5 / f w ), 7r a denotes nine pseudoscalar fields (770,7?, rjs) with 
i =1,...,4 and M(q 2 ) is the constituent mass. This constituent quark mass, which vanishes at 
large momenta and is frozen at low momenta at a value around 300 MeV, can be explicitly 
obtained from the theory but its theoretical behavior can be simulated by parameterizing 
M(q 2 ) = m q F(q 2 ) where m q ~ 300 MeV, and 



F{q 2 



A 2 



_, 1 

_> 



(4) 



A 2 + q 2 

The cut-off A fixes the chiral symmetry breaking scale. 

The Goldstone boson field matrix [7 75 can be expanded in terms of boson fields, 

= 1 + -% 5 A a 7r a - r^TrV + - (5) 

The first term of the expansion generates the constituent quark mass while the second gives 
rise to a one-boson exchange interaction between quarks. The main contribution of the 
third term comes from the two-pion exchange which has been simulated by means of a 
scalar exchange potential. 

In the heavy quark sector chiral symmetry is explicitly broken and this type of inter- 
action does not act. However it constrains the model parameters through the light meson 
phenomenology and provides a natural way to incorporate the pion exchange interaction in 
the DD* dynamics. 



Beyond the chiral symmetry breaking scale one expects the dynamics to be governed by 
QCD perturbative effects. They are taken into account through the one gluon-exchange 
interaction 33[ derived from the lagrangian 



C gqq = iy/4:7ia s ip-ff,G^X c ip , (6) 

where A c are the SU(3) color generators and G£f the gluon field. 

The other QCD nonperturbative effect corresponds to confinement, which prevents from 
having colored hadrons. Such a term can be physically interpreted in a picture in which the 
quark and the antiquark are linked by a one-dimensional color flux-tube. The spontaneous 
creation of light-quark pairs may give rise at same scale to a breakup of the color flux-tube 



34j. This can be translated into a screened potential [35| in such a way that the potential 



saturates at the same interquark distance. 

VcoNifij) = {~a c (1 - e-^) + A}(A^ ■ X /) (7) 

where A is a global constant to fit the origin of energies. At short distances this potential 
presents a linear behavior with an effective confinement strength a = —a c fi c (X c i ■ X c j) while 
it becomes constant at large distances. It has been shown that this form of the potential 
is important to explain the huge degeneracy observed in the high excited light meson spec- 



trum [36| and turns out to be very important for the correct assignment of J PC = 1 
charmonium states js^]. Explicit expressions for all these interactions are given in \^\ . 

All these ingredients are needed to explain the hadronic phenomenology. Apart from the 
obvious confinement potential, gluon exchange is demanded from the hyperfine splitting in 
charmonium. Moreover, pion exchange is one of the best established interaction in nature 
being its parameters constraint by a huge amount of experiments. When Goldstone boson 
exchanges are considered at the quark level together with the OGE the possibility of double 
counting emerges. This problem has been studied in the literature concluding that the pion 
can be safely exchanged together with the gluon [ssl . 

Constituent quark models are also criticized because they only incorporate a limited sector 
of the Fock space. In particular its applicability to high excited states may be questionable 
as more thresholds open up. In our case the parameters of the model has been fixed in 
the low lying part of the spectrum where these effects are more easily incorporated into 



6 



them. Furthermore the main contribution of the open channels are taken into account by 
the screened confinement potential. 



B. The coupled channel approach 



To model the 1 ++ cc system we assume that the hadronic state is 

l*> = E c «l^> + Y.Xp{P)\(t>Kh(t>M 2 P) 

a 



(8) 



where \ip a ) are cc eigenstates of the two body Hamiltonian, (j)^ are cn (en) eigenstates 
describing the D (D) mesons, \4>m 1 4>m 2 (3) is the two meson state with (3 quantum numbers 
coupled to total J PC quantum numbers and Xp{P) ls the relative wave function between the 
two mesons in the molecule. As we always work with eigenstates of the C-parity operator 
we use the usual notation in which DD* is the right combination of DD* and D*D. 

The coupling between the two sectors requires the creation of a light quark pair nn. 
Similar to the strong decay process this coupling should be in principle driven by the same 
interquark hamiltonian which determines the spectrum. However Ackleh et al. 



have 



shown that the quark pair creation 3 P model 



401 ]. gives similar results to the microscopic 



calculation. The model assumes that the pair creation Hamiltonian is 



% — g I d 3 x ip(x)ip(x) 



(9) 



which in the non-relativist ic reduction is equivalent to the transition operator 41] 



T = -3V2iJ2 / d 3 pd 3 p'6 {3) (p + p') 



yi 



p — p 



bi( P )dt(p') 



C=l,j=0,S=l,J=0 



(10) 



where fx (u — Jj) are the quark (antiquark) quantum numbers and 7' = 2 5 / 2 7r 1 / 2 7 with 
7 = ^ is a dimensionless constant that gives the strength of the qq pair creation from the 
vacuum. From this operator we define the transition potential Vp a (P) within the 3 Pq model 



as 



42 1 



(Mm^\T\^) = PV/s a (P) 5 i3) (P cm ) 
where P is the relative momentum of the two meson state. 



(11) 
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Using the wave-function from Eq. ([8]) and the coupling Eq. ( TTTj) we arrive to the coupled 
equations 

M Q c Q + W V a p{P) X p{P)P 2 dP = Ec a 

p 

E / H^ M *(P',P) X p(P)P*dP + j:Ve, a (P')c a = E X p'(P') (12) 

P J a 

where M a are the masses of the bare cc mesons and H^ M2 is the RGM Hamiltonian for 
the two meson states obtained from the qq interaction. 

Solving the coupling with cc states we finally end up with an Schrodinger type equation 
for the relative wave function of the two meson state 

E / (H^ M *(P / ,P)+V p % f (P',P))xp(P)P 2 dP = Exe'(P') (13) 

P 

where 

Vp>p \P iP) - 2-j E _ Ma ^) 
is an effective interaction between the two mesons due to the coupling with intermediate cc 
states. 

In this way we study the influence of the cc states on the dynamics of the two meson 
states. This is a different point of view from the usually found in the literature where the 
influence of two meson states (in general without meson-meson interaction) in the mass 



and width of cc states is studied [42j . Our approach allows to generate new states through 
the meson-meson interaction due to the coupling with cc states and to the underlying qq 
interaction. As we will see the renormalization effects of the cc mass due to this channel is 
small. 

The cc probabilities are given by 

ca = — E / V a p(P) X p(P)P 2 dP (15) 
with the normalization condition 1 = J2a \ c a\ 2 + Y^piXplXp) • 

C. Flatte parametrization 

In order to compare the predictions of our model with the recent Belle and BaBar ex- 
perimental data we obtain from Eq. ( TT2l) a Flatte-like parametrization of the DD* near 
threshold amplitude following Ref. [43j. We remind here the main ideas. 
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From Eq. ( !T2|) . and neglecting the DD* interaction, one can easily derive the DD* 
scattering amplitude 

V 2 (P) 

F^(P,P;E) = -nfiY- ' —r (16) 

where the function ^d* (-^) i s given by 

r V 2 (P) 
g a DD ,(E) = V / -jr 2 ^ P 2 dP. (17) 

For small binding energies e = + Md* — i? it can be expanded as 

g% D .(E) = E a DD « + l -Y DD , + 0(4/A/A 2 ) (18) 

where 

_ poo 

E<* DD * =2^ V$ a {P)dP (19) 
= 2tt^V^(0)P (20) 

and A ^> e is the characteristic scale of the V a p production amplitude which may correspond 
to the scale of the quark wave function and it's assume to be much bigger than the binding 
energy of the physical state. 

A straightforward generalization to include the DD* charged states and other channels 
gives the expression for the near threshold DD* scattering amplitude 

i r DD * 



DD* 



(21) 



2P E - E f + § (T D o D , + r £)+ £)*- + T(E)) + C(Ve/A 2 ) 

where T(E) accounts for the width due to other processes different from the opening of the 
near DD* threshold. Eq. (|2ip corresponds to a Flatte parametrization with 



D{E) = E — Ef + Z -{r DOD *o + T D+D ^ + T(E)) + (9(4/i 2 e/A 2 ). (22) 



Now assuming, as in Ref. [28(|, that the short range dynamics of the weak B — > iOT(3872) 
transition can be absorbed into a coefficient B we are able to write the differential rates in 
the Flatte approximation as 

dBrjB -» KD°D*°) 1 jWg) 

dE ~ ®2tt \D(E)\ 2 • {A6) 
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Mab M- al 

FIG. 1: Diagrams included in the quark rearrangement process DD* — > pJ/^p- 



The analysis of the B — > fCX(3872) — > K-k + -k^ J /ip data is more involved because we 
have to calculate the DD* — > 7T + 7T~J/if) transition amplitude. 

This can consistently be done in our formalism assuming that the process takes place 
through the DD* components of the X(3872) which decays in pJ/ip and then into the final 
7r + 7r~ J /ip states. The decay width of the process is given by 

IVw/* = E t maX dk ~ rr f2" \M J x L _ pJ/4) (k)\ 2 . (24) 

JL Jo {M x -E p - Ej/^f + L -f 

The amplitude -MjcUpj/ip * s calculated in our model by the rearrangement diagrams of 
Fig.[TJ averaged with the DD* component of the X(3872) wave function. The rearrangement 
diagrams are calculated following Ref. [44|. The amplitude is given by 

Mfi = E _Ma (25) 

i=a,a;j=b,b 

where 

Ma{P t P) = {4>M[ 0Afj | 1 0jvfi 4 > M2 ) (£m( I @ij FC I £mlM 2 ) ( 2 6) 

and the orbital part can be written as (e.g. for the case (ij) = {ab)) 
{^m^m^H^mAmi) = J d 3 P M[ d 3 P M >d 3 P Ul d 3 P M2 Pm^Pm'W^PmOKPk-Pmi) 

5{P m > ~ Pm 2 - {P' - P))H{- l -{P Ml + P M2 ) + P M , + 1 -{P' - P)) 
0Mi(-Pmi)0m 2 (-Pa/ 2 )- (27) 
10 



The spin- flavor-color matrix elements are taken from Ref. 44]. 

Once the decay width T n + n -j/^ is calculated, the differential rate is given by 



dBr{B ->■ Kir + n-J/ip) 
dE 



B 



1 ^TT+TT J/lp (E) 



2tt 



(28) 



\D{E)\ 2 ' 

In order to compare with the experimental data we determine the number of events distri 
butions from the differential cross section 



NbJI\E) = 2.5[MeV] 
NEl^\E) = 2.0[MeV] 
1&(E) = 5[MeV] 
NZEZiE) = 2.0[MeV] 



131 \ dBr(B Kn+n-J/ip) 



8.310" 1 
48.3 



dE 

dBr(B -> KD°D°7r ) 



0.73 10~V dE 
93.4 \ dBr(B -> Kn+n~J/i;) 

8.410-V 



33.1 



1.67 10~ 4 



dE 

dBr(B ->■ KD°D*°) 



In all reactions a background is taken into account modelled as in Ref. [28 
KD°D°7T the D°D*° signal interferes with the background and so a phase 



(29) 
(30) 
(31) 
(32) 

For the B ->■ 

^BeZZe = qO &nd 



LBaBar 



324 have been introduced. Also the experimental branching ratio B(D 



*o 



D°tt°) = 0.62 is introduced. We use a value for B 
one used in Ref. 



3.510 which is in the order of the 



28] 



III. RESULTS 



Numerical methods 



To found the quark-antiquark bound states we solve the Schrodinger equation using the 
Gaussian Expansion Method 45j. In this method the radial wave functions solution of the 



Schrodinger equation are expanded in terms of basis functions 

R a (r) = £ 6£0S(r) 



(33) 



n=l 



where a refers to the channel quantum numbers. The coefficients b" and the eigenenergy E 
are determined from the Rayleigh-Ritz variational principle 



E 

n=l 



(34) 



11 



where the operators T% n and N°, n are diagonal and the only operator which mix the different 
channels is the potential Vff^ . 

To solve the four body problem we also use the gaussian expansion of the two body wave 
functions obtained from the solution of the Schrodinger equation. This procedure allows us 
to introduce in variational way possible distortions of the two body wave function within 
the molecule. Using these wave functions Eq. (13) reduces to a matrix equation by Gauss 
integration. 

A crucial problem of the variational methods is how to choose the radial functions (f>ni( r ) 



in order to have a minimal, but enough, number of basis functions. Following |45] we employ 
gaussians trial functions whose ranges are in geometric progression. The geometric progres- 
sion is useful in optimizing the ranges with a small number of free parameters. Moreover the 
distribution of the gaussian ranges in geometric progression is dense at small ranges, which 
is well suited for making the wave function correlate with short range potentials. The fast 
damping of the gaussian tail is not a real problem since we can choose the maximal range 
much longer than the hadronic size. 

B. Results 

The calculation is parameter free since all the parameters are taken from the previous 
calculation {jy], including the 7 = 0.26 parameter in Eq. (flOl) . This value was fitted to 
the reaction ^(3770) — > DD which is the only well established charmonium strong decay. 
This way to determine the value of 7 might overestimate it since the ■0(3770) is very close 



to the DD threshold and FSI effects, which were not included, might be relevant 46]. 

We first perform an isospin symmetric calculation including 3 Si and 3 Di DD* partial 
waves and taking the D and D* masses as average of the experimental values between 
charged states. If we neglect the coupling to cc states we don't get a bound state for the 
DD* molecule in the 1 ++ channel, neither in the 1 = nor in the 1 = 1 channels. The 
interaction coming from OPE is attractive in the 1 = channel but not enough to bind the 
system, even allowing for distortion in the meson states. 

Now we include in the 1 = channel the coupling to cc states. The most relevant are the 
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M (MeV) 

A VA \ A V A. V / 


cc( l 3 Pi ) 


cc(2 3 Pi 1 


D°D*° 


D ± D* T 




3936 


0% 


79% 


10.5 % 


10.5 % 


A 


3865 


1% 


32% 


33.5% 


33.5% 




3467 


95% 


0% 


2.5% 


2.5% 




3937 


0% 


79% 


7% 


14% 


B 


3863 


1% 


30% 


46% 


23% 




3467 


95% 


0% 


2.5% 


2.5% 




3942 


0% 


88% 


4% 


8% 


C 


3871 


0% 


7% 


83% 


10% 




3484 


97% 


0% 


1.5% 


1.5% 



TABLE I: Masses and channel probabilities for the three states in three different calculations. The 
first three states are found when we perform and isospin symmetric calculation with a value of 7 
fit to the decay ^(3770) —> DD. The second three states shows the effect of isospin breaking in 
the DD* masses. The last three states correspond to a value of 7 = 0.19 that fits the experimental 
mass of the A(3872). The probability is shown as zero when it is less than 0.5%. 

1 ++ ground and first excited states with bare masses within the model given by 

cc(l 3 Pi) ->• M = 3503.9 MeV 

(35) 

cc(2 3 P!) -> M = 3947 .4 MeV. 

The results of this calculation are shown in part A of Table HI We find an almost pure 
cc(l 3 Pi) state with mass 3467 Me V which we identify with the Xci(l-P) and two states 
with significant molecular admixture. One of them with mass 3865 MeV is almost a DD* 
molecule bound by the coupling to the cc states. The second one, with mass 3936 MeV", is 
a cc(2 3 Pi) with sizable DD* component. We assign the first state to the X(3872), being 
the second one a candidate to the X(3940). We have also analyzed the effect of higher bare 
cc states finding a negligible effect on the mass and probabilities that will not change the 
above numbers. 

Coexistence of the uJ/ip(I = 0) and pJ/^(I = 1) decay modes strongly suggest a large 
isospin mixing. However the relative branching fraction of both modes can be misleading 
with respect to the absolute magnitude of the isospin mixing in X(3872) due to the phase 
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space suppression of the uJ/ip channel against the pJ/ip one. In fact if we assume that 
X(3872) is a D°D*° molecule, the ratio B{ ^x(S^Zl+l-)^) would be a factor 20 smaller 
than the experiment due to the different phase space. 

It is clear that we need charged components in the wave function but with a different 
weight with respect to the neutral component. This rules out the intuitive idea of the 
dominance of the loosel y b ound neutral component. The clarification of this puzzle has 



been nicely done in Ref 



47|. 



To introduce the isospin breaking in our calculation we turn to the charge basis instead 
of the isospin symmetric basis with the transformation 

{D^D**) = -=(\DD*I= 0) - \DD*I = 1)) (36) 
v2 

\D°D*°) = —= (\DD*I = 0} + \DD*I = 1}) (37) 
a/2 

writing our isospin symmetric interaction on the charged basis. We now explicitly break 
isospin symmetry taking the experimental threshold difference into account in our equations 
and solving for the charged and neutral components. Of course, if we don't break it explicitly 
we recover our previous result as a bound state in the 1 = sector. Now we get again three 
states being the main difference in the DD* molecular component. The masses and channel 
probabilities are shown in part B of Table [B We now get a higher probability for the D°D*° 
component although the isospin component still dominates with a 66 % probability and a 
3% for isospin 1. 

Having in mind that the 3 Po model is probably too naive and we might be overestimating 
the value of 7, we show in Fig. [2] the variation of the X(3872) mass with it. We can see 
that it is possible to get the experimental binding energy with a fine tune of this parameter. 
Using 0.6 MeV as the binding energy we get a value of 7 = 0.19, 25% smaller than the 
original. The results are shown in part C of Table [H Now the D°D*° clearly dominates 
with a 83 % probability giving a 70 % for the isospin component and 23 % for isospin 1. 



Of course, as the isospin breaking is a threshold effect [25j, it grows as we get closer to it 
as can be seen in Fig. [3] where we show the probabilities of the different components for the 
state X(3872). 

In Fig. H] we compare our results with the B — > KD°D°7c° data from Belle (a) and 
B — > KD°D*° data from BaBar (b). The same comparison is done in Fig. [5] for the B — > 
K7r + 7T~ J/\& data from Belle (a) and BaBar (b). In all figures the dashed lines shows the 
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FIG. 2: Mass of the X(3872) as a function of the strength 7 of the 3 Po model. The isospin 
symmetric calculation is shown in figure (a) and the isospin breaking in figure (b). Dotted lines 
show the threshold positions for the DD* average in figure (a) and D°D*° and D ± D* T in (b). 
The solid lines shows the full result and the dashed lines turning off the DD* interaction. 



results without resolution functions. The solid line gives the result using the resolution 
functions as in Ref. 28]. All the resolution functions are those given by Belle 13] and 



BaBar 



19| collaboration with the exception of the BaBar DD* resolution where we use the 



prescription from Ref. 



We find a good description of the Belle B — > KD°D°tt° data whereas the agreement is 
poor in the case of the BaBar data. It is important to notice that in the Belle analysis 
the mass of the X appears as 3872 MeV while in the BaBar data the resonance is located 
3 MeV above. The BaBar mass value does not coincides with the mass of the X obtained 
in our calculation which may be the reason for the disagreement. 

The B — > K7T + 7T~J/ty data are equally well described for the Belle and BaBar experi- 
ments. In this case both Collaborations give similar values for the mass of the resonance, 
namely 3871.4 Me V, which are in much better agreement with our result. 
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FIG. 3: Probability (in %) of different components as a function of the binding energy when we vary 
the 7 parameter of the 3 Pq model. The solid line gives the D°D*° probability, the dashed-dotted 
the D ± D* T , the dashed the cc(2 3 Pi) and the dotted the cc(l 3 Pi). 

IV. SUMMARY. 



As a summary, we have shown that the A" (3872) emerges in a constituent quark model 
calculation as a dynamically generated mixed state of a DD* molecule and Xci(2P). Al- 
though the cc mixture is less than the 10% it is important to bind the molecular state. 
This result is in agreement with the analysis of Ref 29) . The proposed structure allows to 
understand simultaneously the isospin violation showed by the experimental data and the 
radiative decay rates. Furthermore, we have demonstrated that this solution explain the 
new Belle data in the D°D°7r° and 7t + tt~J/^ decay modes and the 7t + 7t~J/^ BaBar data. 
The original Xci(2P) state acquires a significant DD* component and can be identified with 
the A(3940). 
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